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ABSTRACT 



The quantity Yx, the product of the X-ray temperature Tx and gas mass M g , has recently been proposed as a robust low-scatter 
mass indicator for galaxy clusters. Using precise measurements from XMM-Newton data of a sample of 10 relaxed nearby clusters, 
spanning a Y x range of 10 13 -10 15 M G keV, we investigate the M 500 -Y x relation. The M soo - Y x data exhibit a power law relation 
with slope a = 0.548 ± 0.027, close to the self-similar value (3/5) and independent of the mass range considered. However, the 
normalisation is ~ 20% below the prediction from numerical simulations including cooling and galaxy feedback. We discuss two 
effects that could contribute to the normalisation offset: an underestimate of the true mass due to the HE assumption used in X-ray 
mass estimates, and an underestimate of the hot gas mass fraction in the simulations. A comparison of the functional form and scatter 
of the relations between various observables and the mass suggest that Yx may indeed be a better mass proxy than Tx or M s ,5oo- 

Key words. Cosmology: observations, Cosmology: dark matter, Galaxies: cluster: general, (Galaxies) Intergalactic medium, X-rays: 
galaxies: clusters 



1. Introduction 

All theoretical approaches characterise galaxy clusters in terms 
of their mass. Models of structure formation predict the space 
density, distribution and physi cal properties of clu sters as a func- 
tion of mass and redshift (e.g. Bertsch ingerll998l) . However, the 
mass is not easily measured. X-ray estimates from the hydro- 
static equilibrium (HE) equation are valid only for reasonably 
relaxed clusters and require temperature profiles of high statisti- 
cal quality; furthermore, the available precision rapidly degrades 
with redshift. Based on the regularity of the cluster population, 
other X-ray observables such as luminosity, temperature Tx, or 
gas mass M g , have been used as proxies for the mass, e. g. to 
const rain cosmological parameters using cluster surveys 



, I2005I) . Studies of cluster formation physics must also rely on 
mass proxies when considering unbiased (i.e., coverin g a variety 
of dy namical states) or distant cluster samples (e.g. Maughan 
|2007|) . The identification of the best mass proxy, and knowledge 
of its exact relation to the mass, are therefore important. 

The most commonly used mass proxy, Tx, is expected to 
be closely related to the mass via the virial theorem. Significant 
progress on the calibration of the local M-Tx relation for re- 
laxed clusters has recently been made, with ex cellent agree- 
ment now achieved between various observations dArnaud et alj 
2005; VikhlininetaD|2006), and comparison between observa- 
tions and numerical models including cooling and g alaxy feed- 
back showing agreement to the ~ 10% level (e.g . iNagai et aU 
l2007bHArnaud et alJl2005l) . lKravtsov et ail (120® recently pro- 
posed a new mass proxy, Yx = 7jrMg,500> where M gj 5oo is the gas 
mass within /?5oo, the radius corresponding to a density contrast 
of 6 — 500. Yx is related to the thermal energy of the gas and is 
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the X-ray analogue of the integrated SZ Comptonisation param- 
eter, Y$z- The numerical simulations of iRravtsov et al.l showed 
that, as compared to Tx or M gi 5oo> Yx is a better mass proxy, in 
the sense that the intrinsic scatter was lower than for any other 
mass indicator, regardless of cluster dynam ical state (similar to 
previous results for the M-Ysz r elation, e.g.. Ida Silva et al.f2 004; 
iMotl et at] 120051; lNagaill2006l) . Furthermore, its evolution ap- 
pears to be close to the standard self-similar expectation. 

In this Letter, we present the Msqq-Yx relation derived from 
precise XMM-Newton data and compare it to the Msoo-^g,500 
and M^-Tx relations. The A/500-lx relation is discussed with 
respect to previous C handra results and theoretical expectations 
(Nag ai etaill200 7b). Other relations between observables, such 
as the variation of the gas mass fraction / g ,5oo with mass, are also 
investigated in order to shed new light on the scatter and slope 
of the various mass-proxy relations. 



2. The data 

2.1. The sample 

The sample comprises ten nearby morphologically relaxed clus- 
ters in the temperature range [2 - 9] keV. We have previously 
used XMM-Newton data to study the structural and scaling prop- 
erties of the total mass ( Poin tecouteau et aill2005t [Arnaud et al.l 
120051) and of the entropy dPratt et al.ll2006l) ; the T x , M 500 and 
/?soo values derived in these papers are used in the present Letter 
(Table [TJi. T he observations and data r eduction steps are fully 
described in Pointe couteau et ail d2005l) . M500 values were de- 
rived from NFW model fits to mass profiles measured down 
to d b s = 600 - 700, except for the two lowest mass clusters 
(£obs ~ 1400), thus the M500 estimat es involve some data ex- 
trapolation. However, as discussed in lArnaud et alj d2005l) . the 
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Table 1. Physical cluster parameters. M500, Mg ^oo and / gi 5oo are the total mass, gas mass and gas mass fraction respectively, within 
the radius /?5oo, inside which the mean mass density is 500 times the critical density at the cluster redshift. Tx is the spectroscopic 
temperature within [0.15 - 0.75]/?5oo and Yx = M g; 5oo.7x- Values are given for a ACDM cosmology with Q m = 0.3, Qa = 0.7, 
Hq = 70 km s _I Mpc~'. Errors are 1 cr. 



Cluster z Tx(keV) M500 (lO 14 M a ) M g , 500 (10" M G ) Y x (10" M keV) / B , 5 oo 

A 1983 O0442 2.18 ±0.09 1.09^| 0^4^™ 1.39 + ^ 0.058 + ™$ 

MKW9 0.0382 2.43 ± 0.24 0.88 + ° : ^ a49 -0Q5 lA9+ -on °- 055 -oo!2 

A 2717 0.0498 2.56 ± 0.06 1,10«| L02 -ol 2 - 60 -flt5 a093 -oo!l 

A 1991 0.0586 2.71 ±0.07 1.203$ L25 -oi 3.39!™ °- 104 -ooh 

A 2597 0.0852 3.67 ± 0.09 2.22+° : ^ 2.51+™ 9 - 21+ Q38 0,113 QMl 

A 1068 0.1375 4.67 ±0.11 3.87 + ° : ^ 3.77 + ° : |° 17 -6 + ofl °- 097+ ooot 

A 1413 0.1430 6.62 ±0.14 4 - 82 -o4n 7 -55^ 50.0+| j °- 157+ nn!4 

A 478 0.0881 7.05 ±0.12 7.57 + [ : | 9.33+f™ 65.8+| : * 0.123+™ 

PKS 0745-191 0.1028 7.97 ± 0.28 7.27+°^ 10 - 71 -o47 85.3^° 0.147* 

A 2204 0.1523 8.26 ± 0.22 8.39^ 10.55^ 87.2^ : 0-126+jj : °j* 



Table 2. Observed scaling relations. For each observable set 
(B, A), we fitted a power law relation of the form B = C{A/Ao) a , 
with A = 5 keV; 4 x 10 13 M o ; 2 x 10 14 M Q keV for T x , M g , 50 o and 
Yx respectively. crj ogjI and o"i ogi ; are the raw and intrinsic scatter 
about the best fitting relation in the log-log plane . The M ^qq-Tx 
relation is the same as that given in lArnaud et al.l d2005b . 
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M500 estimates rely solely on the physically and observationally- 
motivated assumption that the best fitting NFW model remains 
valid between 5 D \, S and 6 = 500, and not on a less reliable ex- 
trapolation of density and temperature profiles. The temperature 
Tx was derived from a single-temperature fit to the integrated 
spectrum in the [0.1-0.5] ^200 aperture, the inner radius defined 
to exclude the cooling core region and the outer radius chosen 
to ensure a sufficiently precise Tx estimate over the whole mass 
range. This aperture corresponds to [0.15-0.75] ^500, while an 
aperture of [0.15-1] R500 is used for the definition of Tx in nu- 
merica l simulations and in the Chandra analysis (Naga ietall 
2007b). For typical decreasing temperature profiles, these Tx 
values are expected to be slightly smaller by 3-6%Q 

The integrated gas mass depends sensitively on the gas den- 
sity at large radius. To compute M gj 5oo, we re-derived the gas 
density profile from the emissivity corrected surface bright- 
ness profiles using the deproj ection and PSF - decon volution tech- 
nique recently developed bv lCroston et al.l d2006l) . This deriva- 
tion is free of any assumption on profile shape, such as power 



10 13 10 14 10 15 

Yx=M g ,5oo T x (h 70 " 5/2 M keV) 

Fig. 1. The A/500-Fx relation as seen by XMM-Newton from 
a sample of 10 local relaxed clusters. The red solid line is the 
best fitting power law and the shaded orange area corresponds 
to the 1 cr uncertainty. The predicted relation from numerical 
simula tions including cooling and galaxy feedback (Naga i~et~ail 
l2007bl) is over-plotted as a green dot-dashed line (true mass) and 
as a green solid line (mass estimated from mock X-ray observa- 
tions and the HE equation). The dotted and dashed bl ue lines are 
the obse rve d relations derived from Chandra data by Naga i et al.l 
(l2007bl) and lMaughanl d2007l) respectively (see text). 



law behaviour at large radius, a feature common to all analyt- 
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ical fitting models used t hus far (e.g. 



1 The difference is 3% for A1413 dArnaud et all 120051) . a 
cluster for which the temperature profiles measured up to R500 
both with XMM-Newton and Chandra are in excellent agreement 
dPratt & Arnaudl 120021 ; IVikhlinin etail l2005h . In the simulations of 
Nagai et all i2007ah . the [0. 15-0.5] R 50 o temperature is 6% higher than 
that in [0.15-1] Rsoo- A smaller difference is expected for the aperture 
used here. 



Point ecouteau et ail 120051; IVikhlinin et al I l2006h . Furthermore, 
the statistical errors ar e readily estimated fr om a built-in Monte- 
Carlo procedure (see Cros ton et al.ll2006l for details). For the 
present sample, there is excellent agreement between the depro- 
jected density profiles an d the analytical model profiles derived 
in ou r previous work dPointecouteau et all 120051; iPratt et all 
2006). The significan t differences are in t he very central regions 
of some clusters (e.g. Crosto net al.l2006l Fig 12) and for A2597 
at large radii, where the deprojected profile is slightly steeper 
than the model profile. The gas mass estimated with the two 
methods differs by less than 3%, except for A2597 (8% differ- 
ence). For all clusters, except for A1983 and MKW9, the sur- 
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Fig. 2. Correlations between X-ray observables. Solid lines: best fitting power law relations. Dashed lines: standard self-similar 
relation (slope indicated in each figure) normalised to data from the th ree most massive clu sters. Green dotted line in top-left panel: 
best fitting power law Msqq-Tx relation for the hot cluster sub-sample ( Arna udet al.l2005h. Green lines in bottom-middle panel: gas 
mass fraction from the Msoo-Mg^oo relation in the numerical simulations of iNagai et alJ (2007b), using the true mass (dash-dotted 
line) and the HE mass (full line). 



face brightness profiles extend at least up to /?5oo, or very close 
to it, so that extrapolation uncertainty is not an issue. For A1983 
and MKW9, the M g ,5oo estimated from extrapolation in the log- 
log plane are 31%°(A1983) and 67% (MKW9) larger than the 
gas mass measured at d> bs ; as shown below, these points do not 
however have a significant effect on the results. 

The resulting M gi 5no, /s,500 and Fx values are listed in 
Table Q] Errors on M g 5oo include both statistical errors and er- 
rors due to uncertainties in ^500, which are summed quadrati- 
cally, with the latter dominating the error budget. 

2.2. Scaling relations 

For each observable set (B, A) we fitted a power law relation of 
the form h(z)"B - C(A/A ) a , where h(z) is the Hubble constant 
normalised to its present value and n is fixed to the expected scal- 
ing with z. The fit was performed using linear regression in the 
log-log p lane, taking into account the errors on both variables 
(FITEXY: |Press et al.ll 1992b . The pivot point A is chosen so that 
the normalisation and slope are nearly independent parameters. 
For the Msqq-Yx relation for instance, the covariance in log(C) 
and a normalised to the product of their standard errors is 0.042. 
The resulting values are given in Table [2j and the various corre- 
lations are plotted in Fig.Q]and Fig. [2] Table|2]also shows the raw 
and intrinsic scatter about the best fitting relations in the log-log 
plane. The raw scatter was estimated using the vertical distances 
to the regression line, weighted by the error. The intrinsic scat- 
ter was computed from the quadratic difference between the raw 
scatter and the scatter expected from the statistical errors. 



The regression method is strictly valid only if the intrinsic 
scatter is negligible as compared to the statistical scatter; in fact 
they are of the same order (Table [2]). We verified that the re- 
su lts are unchanged using the variation of the method discussed 
in lPratt et alJ (2006). Finally, the Msoo-Ix relation is robust to 
exclusion of A1983 and MKW9, for which data extrapolations 
were required (see above): the differences are at the +0.4 cr and 
-0.2 cr levels for the slope and normalisation, respectively. 



3. Discussion and conclusions 

3. 1 . Comparison with theoretical predictions 
The slope of the observed relation: 



Kz) 2/5 M 50l) = 10 



14.556+0.015 



2x 10 14 M o keV 



0.548±0.027 



h 7 - ( j Mo(l) 



is slightly smaller than the standard self-similar value (or = 3/5), 
at the 1 .9 cr significance level, consistent with the A/500-Mg 300 
and M500— Tx relations (Tableland Fig.|2]top panel). The M500- 
M g 500 relation is shallower than expected, reflecting the increase 
in gas mass fraction with mass (Fig. [2] bottom-middle panel), 
while the M^oo-Tx is steeper. At a given mass the gas mass is 
smaller and the temperature is higher, leading to a partial can- 
cellation in the product Fx = M g Tx (see also below). 

The observed normalisation is ~ 20% smaller over the whole 
Fx range than that derived from num erical simulations includ- 
ing cooling and galaxy feedback (Nagai ^taD l2007bh . while 
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the observed slope is consistent with the predicted slope, a - 
0.568 + 0.006, within the 1 cr error (Fig. [I). Better agreement is 
obtained with the simulated M?3.— Yx relation, where M?3. is the 
mass estimated from mock X-ray observations and the HE equa- 
tion. Although the predicted slope, a = 0.596 + 0.010, is slightly 
higher, the difference in normalisation drop s to ~ 8% ( 2.4 cr) a t 
Y x = 2 x 10 14 M o keV. As discussed by iNagai et ail d2007bl) . 
the offset in normalization, also observed with Chandra data, 
may arise from an underestimate of the true mass by the HE 
equation, perhaps due to residual non-thermal pressure support. 
These numerical simulations also predict a hot gas mass frac- 
tion systematically smaller than observed (Fig.|2]bottom-middle 
panel). The difference is smaller for simulated / g ,5oo using M?£. 
and again could be due, in part, to biases in X-ray mass esti- 
mates. Nevertheless, there may also be an underestimate of / gi 5oo 
in the simulations, possibly due in par t to over-condensa tion of 
hot gas into the cold dense phase (Nag ai et"aill2007bl) . This 
would contribute to the offset, by shifting the M^oo-Yx relation 
to the left in the log-log plane. Finally, as the normalization de- 
pends on r£ 6 , the difference in the exact definition of Tx (see 
Sec. 12.11 ) could contribute by < 4% to the offset. 



that the logarithmic scatter in the M g 5oo-Tx relation is 1/3 of the 
scatter in the M gj 5oo-7x relation for M g 5oo °c T^ (Table |2). 

In terms of observed scatter in the relation with mass, Fx 
thus does not appear to be a better proxy than Tx, and is only 
slightly better than M gl 5oo. However we caution against over- 
inerpretation. Firstly, the present results are for relaxed clusters 
only: with the current data w e cannot check if the s catter is in- 
sensit ive to dynamical state (iRravtsov et al J 120061: iPoole et ail 
|2007|) . Secondly, the scatter estimates should be confirmed us- 
ing larger cluster samples with stricter selection criteria. 

However, in terms of functional dependence with mass, Y x 
is clearly a better proxy than M gl 5oo: it is better fitted by a 
simple power-law, and has a slope closer to the standard self- 
similar value (Table|2]i. Furthermore, although the quality of the 
power law fits to M^oo-Tx and A/500-lx are formally similar 
(X 2 /d.o.f ~ 13/8), with similar (~ 2 cr) deviations from the stan- 
dard slope, there is some indication that Tx is also a better proxy 
than Tx in this regard. The slope of the Ms qq-Tx relation may 
depend on mass range (lArnaud et al.l l2005). reaching the stan- 
dard value when cool clusters are excluded, but the slope of the 
Msoo-Yx relation remains stable in that case (0.7 cr difference). 



3.2. Comparison with Chandra results 

Our M^qq- Yx relat ion is very similar to that derived by 
Nagai et al. (2007b) from the Chandra data presented in 
Vikh linin et alj d2006l see our Fig. []]). The slope a = 0.526 + 
0.038 is consistent with our value, a - 0.548 + 0.027, and 
the normalisation at Y x = 2 x 10 14 M o keV, M 500 = 3.82 X 
I0 u h^ M Q , is higher than our value, (3.60 + 0.13) x 10 14 M G , 
at only the 1 .6 cr level. Even bett er agreement is ob tained with 
the best fitting relation quoted by Maughan (120071 dashed line 
in Fig. []]), derived from the same data excluding the lowest 
mass cluster (A. Vikhlinin, priv. communication). Here the slope 
(a = 0.564) is closer to the self-similar value, as we have found, 
and the difference in normalisation is less than 5% over the 
whole mass range. 

3.3. Comparison of mass proxies for relaxed clusters 

For relaxed clusters. [Kravtsov et alj (|2006) found similar scat- 
ter in the Msoo-^x and M5oo-M g , 5 oo relations (cr log = 0.055 and 
0.047 respectively), but two times less scatter in the M^qq-Yx 
relation (0.022). We can compare with the present data, the sta- 
tistical quality allowing us to estimate the intrinsic scatter for 
the first time. The scatter (Table |2]i is the same for the M500-TX 
and Msoq-Tx relations (cr lo£ , ; = 0.039) and slightly larger for 
the M5oo-M g ,5oo relation (cn og ,i = 0.044). The latter may reflect 
that the Msoo-Mg 500 relation is not actually a power law: the gas 
mass fraction appears constant at Mgoo > 2 - 3 x 1O 14 M , with 
a progressive drop at lower mass (Fig.|2]bottom-middle panel). 

In fact the behaviour of / g ,soo appears to be the primary fac- 
tor driving the scatter in the Msm-Yx relation. The M gi 5oo-Tx 
relation is extremely tight (Fig.|2]and Tabled, being well fitted 
by a power law with no measurable scatter, in spite of the pre- 
cision of the data. Since M500 = A/ gj 5oo/./g,500> me scatter in the 
Msq[)-Yx relation simply reflects the scatter in the / gi 5oo-ix rela- 
tion (cf. top and bottom left panels of Fig. [2j. This scatter could 
arise from true scatter in / g; 5oo and/or scatter in the X-ray mass 
to true mass ratio, e.g., due to variations in the magnitude of 
nonthermal pressure support. Note that a low-scatter correlation 
between M g 5oo and Yx is expected: it is straightforward to show 



3.4. Concluding remarks 

Our results suggest that the various mass scaling relations might 
be better understood by considering the gas thermal energy (Y x ) 
and mass (/ g ,5oo) as its most fundamental properties. Let us sup- 
pose that the thermal energy content of the gas is the quantity 
most closely related to the mass (i.e. the best mass proxy is in- 
deed Tx), and that its relation with mass has a quasi-standard 
slope. Let us further note that the gas mass fraction appears con- 
stant at high mass, with a progressive decrease below a 'break' 
mass (reflecting gas loss or incomplete accretion in low mass 
systems due to non gravitational effects). Since varies as 

(Ml T^ /5 ) 5/2 / g 3/2 , one then expects a steepening of the M 500 -T x 
relation at low mass, with a standard slope at high mass. 

A deeper understanding of the mass scaling relations will 
come from the X-ray study o f larger unbiased samp les of local 
clusters, such as REXCESS (Bohringeret al. 2007), combined 
with lensing data. This is necessary to ascertain the dependence 
of the Msoo-Yx relation on dynamical state, and to calibrate its 
normalisation and slope. This step is essential because the use 
of Yx as a mass proxy, as in the case of Tx, requires a detailed 
understanding of non-gravitational effects, in particular of the 
impact of cooling and feedback on the fraction of primordial gas 
that remains in the gravitationally bound hot phase. Precise mea- 
surements at z — are needed to constrain models, on which one 
must rely for high z studies. Significant progress is also expected 
from forthcoming SZ data (e.g from the Planck Surveyor all sky 
survey), especially if combined with XMM-Newton or Chandra 
data, which will allow a full study of the M-Y^z relation. 

Acknowledgements. We thank A. Kravtsov and D. Nagai for useful comments 
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